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Abstract 

We use Cosmic Microwave Background data from the WMAP, SPT, BICEP, and QUaD exper- 
iments to obtain constraints on the dark matter particle mass m x , and show that the combined 
data requires m x > 7.6 GeV at the 95% confidence level for the XX ~^ bb channel assuming s— wave 
annihilation and a thermal cross section (a a v) = 3 x 10 -26 cm 3 /s. We examine whether the bound 
on m x is sensitive to as measurements made by galaxy cluster observations. The large uncertainty 
in eg and the degeneracy with Q m allow only small improvements in the dark matter mass bound. 
Increasing the number of effective neutrino- like degrees of freedom to 2V e g = 3.85 improves the 
mass bound to m x > 8.6 GeV at 95% confidence, for the XX ~~ ^ bb channel. We also study models 
in which dark matter halos at z < 60 reionize the Universe. We compute the Ostriker-Vishniac 
power resulting from partial reionization at intermediate redshifts 10 < z < 60, but find the effect 
to be small. We discuss the importance of the large angle polarization as a complementary probe 
of dark matter annihilation. By performing Monte Carlo simulations, we show that future experi- 
ments that measure the EE power spectrum from 20 < / < 50 can exclude m x ~ 10 GeV at the 
2 (3) a level provided the error bars are smaller than 4 (3) x cosmic variance. We show that the 
Planck experiment will significantly improve our knowledge of dark matter properties. 
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I. INTRODUCTION 



The nature of dark matter remains one of the biggest puzzles in cosmology. Observational 
evidence from large scale structure observations, gravitational lensing, the cosmic microwave 
background (CMB), and galaxy rotation curves have confirmed the existence of dark matter 
if general relativity is the correct theory of gravitation. A well motivated dark matter 
candidate is the Weakly Interacting Massive Particle (WIMP), a good example of which is 
the Lightest Supersymmetric Particle. The dark matter direct detection experiments DAMA 
[H |2], CoGeNT [3l H], and CRESST [5], have obtained results that are consistent with the 
presence of WIMP dark matter in the Galaxy, and tentatively suggest a low dark matter 
mass ~ 10 GeV. While these results are not without controversy, it is intriguing that WIMPs 
in the Galaxy may have been detected. For WIMP masses m x ~ 10 GeV, one may test 
dark matter properties by indirect detection experiments. In this article, we will restrict 
our discussion to the CMB, and how precise measurements of the CMB power spectra may 
constrain dark matter properties. 

Dark matter annihilation results in energy being released in the form of standard model 
particles, some of which is absorbed by the surrounding gas, causing the gas to heat and 
ionize. The excess free electrons scatter CMB photons coming to us from the surface of 
last scattering. Thomson scattering of CMB photons leads to partial homogenization of 
temperature, resulting in a damping of the CMB power spectra on small scales. Thomson 
scattering also polarizes the CMB, which results in a boost in the EE polarization power 
spectrum on large scales. 

Several authors have studied the impact of particle decay and annihilation on the CMB 
power spectra. Ref. •(]] and [7J studied the effect of decaying dark matter on the ionization 
history of the Universe and the effects of a large optical depth on the CMB. Authors jH] 
showed that CMB polarization could detect or constrain dark matter annihilation. Ref. [9] 
and [10J studied the annihilation of very light dark matter particles at early times. Ref. [11] 
examined the impact of dark matter annihilation including the clumping of dark matter, 
while [12] performed a more detailed analysis of dark matter annihilation in halos com- 
posed of light dark matter particles, and the resulting increase in optical depth. Further 
constraints on dark matter properties based on partial reionization and optical depth calcu- 
lations were obtained by [13] , and [H] , while [15] and [16] examined leptonically-annihilating 
dark matter models. Authors [IT] and [18] studied how future large angle CMB polariza- 
tion measurements could detect light dark matter, and ways in which reionization from dark 
matter annihilation could be distinguished from reionization from baryonic sources. Authors 
[19] performed detailed computations of energy absorption at high redshifts. Dark matter 
constraints based on CMB measurements were obtained by [T9H22] . More recently, analysis 
of WMAP data by [23] and WMAP+ACT data by [21] have placed stringent constraints on 
WIMP dark matter with mass m x < 10 GeV. Authors [22] performed a principal component 
analysis of energy injection due to dark matter annihilation, with relevance to current and 
future CMB observations. 

In this paper, we build upon the existing literature in many ways. In Section II, we 
provide an introduction to WIMP annihilation and how it influences the reionization history 
of the Universe, which in turn influences the CMB anisotropies. In Section III, we perform a 
maximum likelihood analysis using the CMB TT power spectrum data from the WMAP [26] 
and SPT [27] experiments, TE data from WMAP, and EE data from the BICEP [2H] and 
QUaD [2H] experiments. We confirm earlier results, and show that WIMP masses m x < 7.6 
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GeV may be excluded at the 95% confidence level by the combined data, for the annihilation 
channel XX bb. Since direct detection experiments favor WIMP masses in this range, it is 
important to investigate whether these bounds may be improved by combining CMB data 
with other experiments. The bound on the WIMP mass is weak because the damping of 
the power spectra (~ exp — 2r) may be largely offset by an increase in A s (the amplitude 
of the primordial curvature power spectrum), or equivalently by increasing the amplitude 
of matter fluctuations (erg). We therefore ask whether the bounds on o~g obtained from 
galaxy cluster observations may strengthen the CMB bound on the dark matter mass m x . 
Unfortunately, the large <jg Q m degeneracy of current galaxy cluster observations does not 
permit a significant improvement on the dark matter mass bound. We show that varying 
the number of effective number of neutrino-like degrees of freedom does have an effect on 
the likelihood function. For the best fit value iV e g = 3.85 obtained by the SPT collaboration 
[27] , we find that the mass bound improves to m x > 8.6 GeV, compared to the bound 
m x > 7.6 GeV for the standard case of N e s = 3.04. 

In Section IV, we revisit partial ionization by dark matter annihilation at intermediate 
redshifts, but this time accounting for dark matter halos. We show that it is possible to 
substantially reionize the Universe with suitable halo parameters. We compute the Ostriker- 
Vishniac / linear kinetic Sunyaev-Zeldovich power spectrum due to partial ionization by 
dark matter annihilation at intermediate redshifts 10 < z < 60. The effect is however 
quite small, and mostly indistinguishable from the standard scenario with no reionization 
at these redshifts. In Section V, we show that the EE power spectrum provides information 
complementary to what may be obtained with the TT power spectrum. We perform Monte 
Carlo simulations, and show that light (~ 10 GeV) dark matter may be strongly constrained 
by measuring the EE power spectrum at multipoles I ^ 20. 



II. DARK MATTER ANNIHILATION 

Let us consider the case where all of the dark matter is composed of Weakly Interacting 
Massive Particles (WIMPs). The lightest WIMP needs to be stable or very long lived in 
order to be a good dark matter candidate. However, dark matter particles can annihilate 
in pairs, producing standard model particles. The dark matter relic density observed today 
Q x h 2 ~ 0.1 is obtained if the annihilation cross section (a a v) ~ 1 picobarnxc, a value often 
considered to be a natural cross section for a weakly interacting particle. We will assume 
here that this is the case. We will also assume that the annihilation cross section is s-wave 
dominated, i.e. (cr a t>) is independent of temperature (or relative velocity v). 

Let p x {z) be the dark matter density at redshift z. The number of dark matter particles 
of mass m x per unit volume is n x = p x /m x . The probability of WIMP annihilation per 
unit time is (a a v)n x , and the number of dark matter pairs per unit volume is n x /2. The 
energy released per WIMP annihilation = 2m x c 2 , and a fraction f em of this energy consists of 
electromagnetic particles which interact with the surrounding gas. Thus, the useful energy 
released per unit volume and per unit time is given by: 

dE Pl(z)(a a v) 

Jem • \*-J 



dtdV ' m x 



Eq. Q assumes that particles and antiparticles are identical. If particles are distinct from 
antiparticles, one must multiply Eq. Q by 1/2. For free dark matter particles (i.e. particles 
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not bound in halos), p x is given by 

pl cc (z) = (i + z) (i p 2 c , lt n 2 x , (2) 

where p cr ; t is the critical density, and Q x is the dark matter fraction today. WIMP annihi- 
lation typically produces a number of particles, each with its own energy spectrum. Let us 
direct our attention to low mass WIMPs, i.e. m x ~ 10 GeV, in agreement with the results 
obtained by the DAMA, CoGeNT, and CRESST direct detection experiments. Authors 
[30H33] have considered light neutralinos within the context of an effective MSSM theory 
that does not assume gaugino mass unification at the GUT scale. For light non-relativistic 
WIMPs, annihilation to W, Z, or Higgs bosons is forbidden by conservation of energy. Direct 
annihilation to photons or gluons is one-loop suppressed if WIMPs have no electric or color 
charge. For the special case of Majorana type fermion dark matter, annihilation to light 
fermions of mass m/ is suppressed by a factor (m//m x ) 2 due to helicity conservation. Thus, 
for Majorana WIMPs, particle annihilation to the most massive fermion pair is favored. For 
WIMP masses m x > 5 GeV, this channel is XX ~^ bb. The dominant leptonic channel is 
XX ~^ t + t~ and the sub-dominant hadronic channel is XX ~^ 

Fig. 1 shows the energy spectrum of neutrinos (summed over flavors) from WIMP an- 
nihilation to bb, cc, and t + t~ from [Ml 135] . The fraction of energy released in the form of 
neutrinos is calculated as: ^ ^ 

U = o / dxx—, (3) 

m x c A J ax 

where x = E/m x is the dimensionless energy, and (dN/dx) Ax is the number of neutrinos 
with dimensionless energy in the range [x, x + Ax]. The spectral function dN/dx is normal- 
ized such that f v < 1. The electromagnetic fraction f em defined in Eq. ([I]) = 1 — f v . For the 
3 channels considered in Fig. 1, we find / em = 0.74 for XX ~~ bb, / em = 0.78 for XX ~ >* c c> 
and /em = 0.65 for XX ~* t + t~ . As an extreme example, we note that WIMP annihilation 
directly to e + e~ would result in / em w 1, however, this channel is unimportant for Majorana 
fermions. Dirac fermion dark matter particles do not suffer from helicity suppression, and 
may annihilate directly to e + e~, but an extra factor of 1/2 (since particles annihilate only 
with antiparticles) in Eq. ([T| means that the bound is weaker by a factor of 2. If there is 
significant asymmetry between particles and antiparticles, the bound is weaker still. If the 
dark matter is the lightest Kaluza-Klein particle (in theories with extra dimensions), the 
annihilation to light fermions is not helicity suppressed. 

Let us first consider dark matter annihilation at high redshifts. The mean free path for 
photons/charged particles scattering with background atoms = A ~ [rib(z)a(E)]~ 1 , where 
nb is the baryon number density at redshift z, a is the scattering cross section, and E is 
the energy of the photon/charged particle (we are assuming here that the wavelength of 
the radiation is much smaller than the Bohr radius, and hence neutral Hydrogen atoms 
may be treated as simply protons and electrons). For radiation to be effectively absorbed, 
it is necessary that A is smaller than the Hubble horizon, i.e A < c/H(z), where H{z) 
is the Hubble parameter at redshift z. For cr(E) ~ the Thomson cross section o"t, we find 
A < c/H for z > 40. Note that prompt photons from WIMP annihilation may have very high 
energies, and hence small cross sections, however, inverse Compton scattering with CMB 
photons at high redshifts is very efficient in producing a flux of lower energy photons. Thus, 
for redshifts z ^> 40, we may make the "on the spot" approximation, i.e. the assumption 
that the energy released by WIMP annihilation at a redshift z is absorbed at approximately 
the same redshift. For redshifts z % 40, we must study the radiative transfer carefully to 
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FIG. 1: Energy spectrum of neutrinos from WIMP annihilation to bb, cc, and t + t~ channels 
(the spectra from v e , and v T channels are added together), x = E/m x , and / em = 1 - fv where 
/„ is the neutrino fraction. For the bb, cc, and t + t~ channels, we find / em « 0.74, 0.78, and 0.65 
respectively. 

obtain the energy absorbed at a given redshift. We will discuss this case later. We note that 
even at high redshifts, only a fraction / a b s ~ 0.85 of the energy is absorbed [19J, where the 
bar indicates an average over redshifts. 

At high redshifts, the energy absorbed per atom per unit time, at a redshift z is given 

by: 



Z(z) = /abs/e 



m x n b {z) 

/abs /cm fh (cr a v) /p C rit\ (^x^) 2 /-, , n3 

m x \h 2 ) n h h 2 [ z) 

3.» x io- £ m r . (4) 



for tt h h 2 = 0.0226, Vt x h 2 = 0.111, (<r A v) = 3 x 10" 26 cm 3 /s, and assuming / abs = 0.85, / em = 
0.74. m is the mean mass per atom rj 1.22 m p where m p is the proton mass and assuming 
24% Helium, h = H / (100 km/s/Mpc) is the dimensionless Hubble parameter at the present 
epoch, Q x and J\ are the dark matter and baryon densities today. 

Energy released by WIMP annihilation partially ionizes gas in the surrounding medium. 
The fraction of the absorbed energy that goes into ionization has been studied by [36-39J. 
Detailed computations [36J indicate that a fraction ?7i on w 0.4 of the absorbed energy is used 
for ionization, for small ionization fractions x- lon . A fraction ^heat ~ 0.2 of the energy results 
in heating, for small x lon , with the remaining energy going into collisional excitations. The 
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ionization fraction x- lon (z) and gas temperature T gauS (z) at a given redshift are computed by 
solving together, the two equations: 

- (1 + z)H(z)^j— - = fi [1 - x ion (z)} i] ion (z)^(z) - n(z)x 2 on (z)a(z) 

-(l + z)H(z)*& = -2T(z)H(z) + 2Vh T {Z) Z(z) + XUZ) [T ^\ ~ T{Z)] . (5) 
dz 3/cb t c (z) 

fi ~ 0.07 eV _1 is the inverse of the average ionization energy per atom, assuming 76% H 
and 24% He, neglecting double ionization of Helium [12]. a is the case-B recombination 
coefficient, T 7 is the CMB temperature, and t c is the Compton cooling time scale ~ 1.44 
Myr [30/(1 + z)] 4 . The last term in the temperature evolution equation accounts for the 
transfer of energy between free electrons and the CMB by compton scattering [40H42] . We 
used x lQn <C 1 and ignored the Helium number fraction in the temperature coupling term. 
In practice, we compute x ion and T gas using a modified version of the publicly available 
RECFAST program HD 12]. 



III. CMB DATA AND CONSTRAINTS ON THE WIMP MASS 

Let us first consider the impact of dark matter annihilation at high redshifts on the CMB. 
At redshifts z > 60, the contribution from dark matter halos is negligible, so we consider 
just the smooth component. Energy from dark matter annihilation partially ionizes the 
surrounding gas. This leads to an excess of free electrons which scatter off CMB photons, 
resulting in partial homogenization of the CMB temperature. The power spectra are thus 
damped relative to the standard ACDM model by a factor exp(— 2rd m ) where Td m is the 
excess optical depth due to free electron-CMB photon scattering. Fig. |2](a) shows the 
evolution of the ionized fraction for the standard scenario without WIMP annihilation (solid 
black curve), as well as models with light WIMPs. At very high redshifts z > 800, particle 
annihilation does not result in excess ionization because the available ionization energy due 
to particle annihilation is small compared to the radiation energy. At lower redshifts, the 
effect of particle annihilation is clear. Fig. [2](b) shows the predicted TT power spectrum 
for the model with no dark matter annihilation, as well as for the light dark matter models. 
The total optical depth t(z) up to a redshift z is defined as: 

f z dz 1 

r(z)=J Q {1 + zf)H{zf f^n h (z f )xUz'). (6) 

At high redshifts, H(z) oc (1 + z) 3 ^ 2 , and hence the integrand oc (1 + z) x l 2 x- mn (z). We see 
that even small changes in Xi on at high redshifts can alter the optical depth r, resulting in 
a damped power spectrum. 

A major difficulty in observing the predicted damping is that the power spectrum oc 
A s exp(— 2r), where A s is the amplitude of the gauge invariant primordial scalar curvature 
power spectrum at a pivot scale fc p i vo t = 0.05 Mpc -1 . Unless the value of A s is fixed by 
other experiments, this degeneracy implies that a damped power spectrum may be partially 
offset by a proportionally larger value of A s , or equivalently, by a larger ag. The degeneracy 
is not exact however, as causality requires that scales larger than the horizon at the epoch 
of particle annihilation be unaffected. Since the horizon ~ (1 + z)~ 3 ^ 2 , particle annihilation 
at high redshifts affects the CMB on moderate and small scales, but leaves the large angle 
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FIG. 2: The effect of WIMP annihilation on the ionized fraction and the CMB TT power spectrum, 
(a) shows the evolution of the ionized fraction with redshift, for the standard scenario without dark 
matter annihilation, as well as 2 models with light WIMPs. (b) shows the TT power spectrum 
for the 3 cases. A larger ionized fraction results in more scattering of CMB photons with free 
electrons, and hence lower power. The exp(— 2r) damping may be partially offset by a larger value 
of A s . 

(small I) power spectrum unaffected. We will use this scale dependence to constrain the 
WIMP mass. 

Fig. [3] shows the data set that we will be considering, (a) shows the TT data from the 
publicly available 7-year data release from the WMAP mission combined with data from 
the SPT experiment. For / > 750, we find that the SPT data has smaller error bars than 
WMAP, while for / < 750, the WMAP data is better. We do not consider data for I > 1500 
due to complications from foreground sources, and secondary processes such as lensing. (b) 
shows the TE power spectrum data from WMAP. (c) shows the EE power spectrum from 
the BICEP (for I < 350) and QUaD (for I > 350) experiments. We use 53 data points from 
the TT power spectrum data set, 33 points from the TE data set, and 15 from the EE data 
set. Also shown in red are the predictions for the standard ACDM cosmology without the 
effect of dark matter annihilation. 

As mentioned previously, we consider the simple scenario in which all of the dark matter 
is made up of WIMPs, and the annihilation cross section (a a v) is given by the thermal value, 
and is assumed to be independent of temperature (i.e. the annihilation is s-wave dominated). 
We obtain constraints on the WIMP mass by performing a maximum likelihood analysis 
using the publicly available CMB Boltzmann code CLASS [43, 44J, varying the following 
parameters: {m x , A s , n s , h, Q^h 2 , Q m h 2 }. n s is the scalar spectral index, h is the Hubble 
parameter today in units of 100 km/s/Mpc. f2b and Q m are the baryon and matter density 
fractions at the present epoch. Single step reionization at = 10.5 is assumed. For 
simplicity, we have only varied the cosmological parameters that are most affected by dark 
matter annihilation. As discussed earlier, the parameter A s is largely degenerate with m x . 
The parameters n s , Q^h 2 and Q m h 2 modify the location and height of the peaks and hence 
may compensate for the effect of m x . Varying h affects the best fit value of A s , as well as 
the value of the likelihood function. VL^h 2 and Q m h 2 also affect the ionization directly (Eq. 
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FIG. 3: The CMB data, (a) shows the TT power spectrum. For I < 750 we use the WMAP 7 year 
data release, while for I > 750, we use data from the SPT experiment, (b) shows the TE power 
spectrum from the WMAP 7 year data set. (c) shows the EE power spectrum from the BICEP 
(for I < 350) and QUaD (for I > 350) experiments. Shown in red is the prediction of the standard 
ACDM model without WIMP annihilation included. Following the WMAP convention, we plot 
(I + 1)Ci/2tt for the TE power spectrum and 1(1 + l)Ci/2ir for the TT and EE power spectra. 



(|4j)). We set the equation of state of dark energy u — — 1. We also assume zero curvature, 
and no running of the spectral index. The primordial Helium fraction is set to 0.24. We 
first assume the number of (massless) neutrino-like degrees of freedom N e g = 3.04. We later 
examine the effect of a larger value of N e g. 

We begin by obtaining the best fit values of {A s , n s , h, fl^h 2 , fl m h 2 } by setting m x = oo 
and minimizing the value of x 2 '- 

^ 2 _ <y ( Cf(data) - (theory) \ 2 

where bC\ is the 1 — a error bar in C\. With the combined TT + TE + EE data, we obtain a 
best fit value Xmm = 92.5/96 d.o.f. (101 data points, 5 fitting parameters with m x fixed to 
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oo), with best fit values {10 9 A S = 2.24, n s = 0.97, h = 0.69, VL m h 2 = 0.1395, VL h h 2 = 0.0225}. 

75.1/77 d.o.f., while for the TT 



With the TT + TE data, we obtain a best fit value x 
data alone, we find x mm = 46.8/48 d.o.f. 
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FIG. 4: The likelihood function and WIMP mass exclusion, (a) shows the normalized likelihood as 
a function of 10 9 A S marginalized over other cosmological parameters, for different WIMP masses. 
As the WIMP mass is increased, the likelihood function peaks for larger values of A s . (b) shows 
the likelihood marginalized over A s , as a function of l/m x . WIMP masses below 7.6 GeV are 
excluded at the 95% level for the bb channel. 

The likelihood function C is defined as — 21n£ = \ 2 + constant. We marginalize the 
likelihood by integrating over variables: 

C(a) = J dbC(a,b). (8) 

Fig. |4]^a) shows the normalized likelihood function for h = 0.69 and marginalized over other 
cosmological parameters, as a function of 10 9 /l s . Shown are likelihood curves for different 
WIMP masses. Small WIMP masses result in a significant damping of the power spectrum, 
and hence require a larger A s to compensate, resulting in curves peaking at larger values 
of A s . The area under each likelihood curve is a measure of how well that model fits the 
data. Fig. |4^b) shows the likelihood function marginalized over {A s , n s , h, VL^li 2 , £l m h 2 }, 
as a function of l/m x . At the 95% confidence level, we are able to exclude a WIMP 
mass m x < 7.6 GeV for the specific channel XX ~^ bb, assuming a thermal annihilation 
cross section and s-wave annihilation, and assuming no prior knowledge of cosmological 
parameters. WIMP exclusion limits for other annihilation channels may also be obtained 
from the above result, since £ oc /em/^x ( see -^1- @)- Thus, for the XX ~ * t~ + t~ channel, 
we have / em = 0.65, which excludes WIMP masses m x < 6.7 GeV. For the extreme case of 
/em ~ 1 obtained for XX ~^ e + e~, we exclude WIMP masses m x < 10.3 GeV. We note that 
ignoring the EE data resulted in slightly stronger bounds on the WIMP mass, and emphasize 
that one must include all the available data when computing dark matter constraints. 
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We have seen in the previous section that WIMP annihilation leads to a damping in the 
small scale TT and TE power spectra, requiring a larger value of A s to compensate for 
the damping. Since A s oc a\ for fixed fl m h 2 , Q^h 2 , and n s , this implies a proportionally 
large value of a 8 , where <r 8 is the RMS matter fluctuation averaged over a scale of 8 Mpc/h. 
A large matter power leads to an abundance of galaxy clusters which may be observed by 
means of X-ray telescopes and by CMB experiments through the thermal Sunyaev-Zeldovich 
effect H5J. 

The number of galaxy clusters N(M\, M2) between masses Mi and M2 is given by the 
formula 

rM 2 rz maK {M) rIM dV 

z min is the lowest redshift probed by the survey, and z msix (M) is determined by the minimum 
flux (or integrated Compton-Y parameter) required in order to make a detection of the 
cluster. dN / dMdV is the number of clusters (or halos) per unit mass per unit volume and is 
obtained numerically [4*6H50] . This function is sensitive to cosmological parameters. dV/dz 
is the cosmological volume-redshift relation. Matching the predicted number of clusters with 
the observations yields a bound on a% which is largely degenerate with f2 m . If cr 8 is accurately 
determined, one may use this bound to limit the variation of A s and hence obtain stronger 
bounds on m x using CMB observations. We convert A s to a 8 using the approximate fitting 
formula given by |51j : 

/Tom; / n^y". /turf" {3A23h) ^ [M™ * (10) 

J V 3.125 V0.024/ V - 14 / V - 72 / °- 76 

G ~ 0.76 is the growth factor today, and / is a "fudge factor". Using the cosmological 
parameters determined by [27J, one obtains agreement between <r 8 and A s , for / = 0.96. 

Fig. [5] shows the Q m as contours for m x = 00 and m x = 10 GeV, marginalized over 
A s ,Q, m h , and VL^h 2 . The solid (red) contours are plotted for h = 0.72, while the dotted 
(blue) contours are plotted for h = 0.69. Also shown are the constraints obtained by 
[52] 153] using Chandra observations of X-ray galaxy clusters (the galaxy cluster contours are 
marginalized over h using a Gaussian prior h = 0.72 ± 0.08). We note that other authors, 
e.g. [SI] find smaller as and larger Q m values by accounting for dynamical dark energy. 
From Fig. [5j we note that the uncertainty in as from galaxy cluster observations is quite 
large due to the degeneracy with Q m . Since this uncertainty is typically larger than the 
uncertainty in as due to dark matter annihilation, only modest improvements in the dark 
matter mass bound are possible with current cluster constraints. 

Up to now, we have assumed N e g = 3.04 which is predicted by the standard model of 
particle physics. Surprisingly, both the SPT [27] and ACT [55] experiments have reported 
deviations from the standard value, using precision measurements of CMB data on small 
scales. The SPT collaboration has obtained a value of N e s = 3.85 ± 0.62, while the ACT 
collaboration finds N e g = 5.3 ± 1.3. Fig. [6] shows the likelihood function marginalized over 
cosmological parameters for the case N eS = 3.85, compared to the standard case iV eff = 3.04. 
The likelihood function is slightly larger at m x = 00 implying a better fit to the data, in 
agreement with the SPT results. The likelihood function also falls off more rapidly, resulting 
in a stronger mass bound m x > 8.6 GeV compared to the bound m x > 7.6 GeV obtained for 
N e fi = 3.04. Ref. [56H58J have pointed out that the late time decay of a particle may mimic 
a neutrino-like degree of freedom. Axion dark matter models also prefer a larger value of 
N eS [591 [60]. 
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FIG. 5: TO erg contours from CMB data, along with the allowed parameter space from galaxy 
cluster observations by Vikhlinin et al. (region within the solid black lines). The solid (red) 
contours assume h = 0.72, while the dotted (blue) contours are for h = 0.69. The cluster results 
from Vikhlinin et al. are marginalized over h assuming a Gaussian prior h = 0.72 ± 0.08. 

' iVeff'= 3.04 — 




FIG. 6: The likelihood function marginalized over cosmological parameters, for iV e ff = 3.04 (solid, 
black) and N e g = 3.85 (dotted, blue). The iV e ff = 3.85 case is sightly preferred over the standard 
model, and falls off more rapidly with decrease in m x . The mass bound for the N e g = 3.85 scenario 
is m x > 8.6 GeV compared to the bound m x > 7.6 GeV for the standard N c q = 3.04 scenario. 
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IV. DARK MATTER HALOS AND SECONDARY CMB ANISOTROPIES 



Up to now, we have ignored the presence of dark matter halos. For weakly interacting 
dark matter, the earliest structures form at a redshift z ~ 60, and have masses > 10~ 6 M Q 
[6T| E2] j although smaller minimum masses may be possible [63J . The minimum halo mass 
is set by the free streaming scale [6TL 162] which in turn depends on the WIMP interaction 
cross section. 



A. Dark matter annihilation in halos 



Let us now look at the effect of dark matter halos on the evolution of the ionization 
fraction. The energy released per unit time and per unit volume dE/dtdV is given by Eq. 
([!]), but to account for halos, we must replace p 2 in Eq. ([2j by [12] : 



p 2 X Mio( z ) = ( X + z f 



dM dnha ^(M,z) 
dM K ' 1 



''200 



drAnr p h {r 



(M,z). 



(11) 



^haio is the comoving number density of halos. The volume integral over the halo is a function 
of both halo mass M and redshift z. We have assumed that halos are truncated at r 2 oo> 
the radius at which the mean density enclosed equals 200 times the background density at 
the time of halo formation. Ph( r ) is the halo density at a distance r from the halo center. 



Eq. ( 11 ) ignores halo-halo interactions, as well as interactions between free dark matter and 



dark matter in halos. 

For dark matter halos with a generalized Navarro-Frenk- White (NFW) [M] profile, we 
have 

PhW = , ; _ PS . ; lfl , (12) 



(r/r s ) a [l + r/r/ : 
and the volume integral over the halo takes the form [12] 



/ 



dr Anr 2 p^(r) 
f(c) 



M p /tt dm 



3 V ^~^m 

c 3 H dx x 2 ' 



m 

2a (l + x) 



-2/3 



[Jodxx 2 



[l + x) 



(13) 



M is the mass of the halo, and c is the concentration parameter, p is the mean density of the 
halo which we set equal to 200 times the cosmological average matter density at the time 
of formation of the halo, e is a dimensionless cutoff scale, required to make the luminosity 
finite for a > 1.5. For the standard NFW profile (a = 1, (5 = 2, e = 0), the parameter /(c) 
is quite small. For example, for c = 10, we find /(10) ~ 150. The value of f(c) rises steeply 
when halos are more cuspy which could occur for e.g., due to adiabatic contraction. For 
a = 1.3, we find /(10) « 10 3 , and for a = 1.5, e = lO" 5 , we find /(10) « 10 4 . 

At moderately high redshifts z > 25, inverse Compton scattering with CMB photons is 
the main mechanism by which high energy charged particles from WIMP annihilation lose 
energy to the CMB, producing a large number of medium energy (< MeV) photons which 
ionize the gas. Here, we will assume that inverse Compton scattering is efficient, and the 
up-scattered CMB photons interact with gas atoms with a cross section ~ <tt- If the cross 
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section may be thought of as independent of energy, the energy absorbed per gas atom, per 
unit time £(z) is given by [18J (compare with Eq. Q): 



CO"t 



-dz' 



l + z 



[l+z')H(z') \l+z> 



( J^L) ( z >\ e ~r(z',z) 

\dtdV) {} 



(14) 



The first term in the integrand comes from the relation dz = —dt H(z)(l + z), while the next 
term accounts for the expansion of the Universe in the time it takes the ionizing radiation to 
reach redshift z having been emitted at redshift z' . t(z', z) is the optical depth from z' to z. 
It can be shown that Eq. (14) reduces to Eq. Q in the limit z' th z, i.e. in the tight coupling 
limit, at high redshifts [18J. Fig. [7| shows the evolution of the ionized fraction Xi on (z), for 
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FIG. 7: Evolution of 

X'mn^z) in the presence of dark matter halos, for different values of /(c). The 
solid (black) curve denotes the standard cosmological recombination model with no dark matter 
annihilation, while the dotted curves include dark matter annihilation. The green curve is plotted 
for free particle annihilation only, while the blue, cyan, and magenta curves account for dark matter 
halos. We set M min = 1(T 9 M Q , and use the Sheth-Tormen mass function. 

different cases. The solid black curve shows the case of no dark matter annihilation, while 
the green curve includes dark matter annihilation (m x = 10 GeV), but ignores halos. The 
blue, cyan, and magenta curves include the effect of dark matter annihilation in halos. In 
all cases, one step reionization at « 10.5 is assumed. For the rather extreme case of 
/(c) = 10 5 , the Universe is nearly completely ionized by z = 10 by dark matter halos alone. 
Early ionization by halos may be tested by optical depth computations [T2tfl4] . However, 
the WMAP constraints on the optical depth are largely dependent on the quality of the 
polarization power spectrum data which has a high signal-to-noise ratio only for small I. It 
is thus possible for very early ionization to be undetected by current data, although some 
constraints may be obtained from TT and TE data jHS]- For small values of /(c), the 
primary CMB is not significantly affected by the inclusion of halos. 
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B. The Ostriker-Vishniac effect 



The Ostriker-Vishniac (OV) / linear kinetic Sunyaev-Zeldovich effect is a second order 
effect caused by the scattering of CMB photons by electrons that have a peculiar motion. 
This effect was first studied by Zeldovich and Sunyaev [15] and by Sunyaev and Zeldovich 
[66] . It was studied in the context of large structure formation by Ostriker and Vishniac 
[67] . and by Vishniac [68]. Here, we merely provide the formulae, following the work of Jaffe 
and Kamionkowski [69] (see also [70]). 

The OV effect has been extensively studied in recent years as a probe of reionization. Par- 
tial reionization due to dark matter annihilation in halos differs from standard reionization 
in 2 aspects: (i) It is gradual and occurs at higher redshifts and (ii) it is more uniform as the 
mean free path of high energy photons/charged particles resulting from WIMP annihiation 
may be larger than the mean separation between halos. To our knowledge, the OV effect 
has not been applied to constrain partial reionization due to dark matter annihilation. We 
now examine whether the OV effect may be used to constrain dark matter annihilation in 
halos. 

The fractional temperature change induced by electrons with a bulk motion along the 
line of sight is [68] : 

AT r 

dt (h ■ v) n e a^e r , (15) 



T 

where h is a unit vector denoting the line of sight, v is the peculiar velocity of the electrons, 
n e is the number density of free electrons, ctt is the Thomson cross section, and r is the 
optical depth. In linear theory applicable at high redshifts, the peculiar velocity v may be 
simply expressed in terms of the matter overdensity. The power spectrum of temperature 
fluctuations will then be expressed in terms of the matter power spectrum, which is assumed 
to take the form given by Eisenstein and Hu pTj. The Ostriker-Vishniac power spectrum 
may then be written as: 

c?r = j^«„ dx G2(z)E2{z)v2(z)s (16) 

G, E, V, and S are dimensionless quantities, x = r] H — T](z)H and 

v(z)Ho = / ^ (17) 



/o E(z') 

and E(z) = \Jd m (l + z) 3 + G and T> are given by: 

G{z) = E{z)^-e- T{z) 

D(z) dD(z) , , 

v ^=dW)^ (18) 

D{z) is the growth function at redshift z. The function S(k) takes the form 

S(k) = -±-( Ho -) 5 S vish (k), (19) 



I67H V c 
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FIG. 8: The linear Ostriker-Vishniac power spectrum, for the different ionization histories shown 
in Fig. [7| Only the rather extreme case of /(c) = 10 5 is distinguishable from the standard 
recombination history without dark matter annihilation. 



where SVish(fc) is the Vishniac power spectrum [681 169]: 

SvishW = k fdy £ dx P{ky)P{k^l + y* - 2xy) ^^Lz^! , (20 ) 



Let us apply Eq. (16) to dark matter reionization. Fig. |8| shows the OV power spectrum 
for each of the dark matter models plotted in Fig. [7j Only the case of f(c) = 10 5 is signifi- 
cantly different. Even in this case, the peak value differs from the standard recombination 
theory by only ~ 20%. One may hope to detect small changes in the OV power spectrum 
at frequencies that minimize the thermal SZ contribution. Even then, the contribution due 
to infrared sources needs to be modeled to high precision. The OV effect may be larger if 
the dark matter halos contain ionized baryons. In that case, one must take into account 
the bulk motion of the clusters themselves. The scenario is then similar to patchy reion- 
ization, and a full non-linear treatment is required. In the linear regime, we do not think 
that the OV effect is particularly useful in constraining dark matter annihilation due to 
halos, except perhaps for the most extreme cases. We however expect these extreme cases 
to be inconsistent with the primary CMB itself, thus limiting the usefulness of computing 
secondary anisotropies. Uncertainties in the reionization redshift also add to the difficulty 
in distinguishing different dark matter models. 



V. THE LARGE ANGLE POLARIZATION POWER SPECTRUM 

As discussed earlier, Thomson scattering of CMB photons by free electrons partially ho- 
mogenizes the CMB temperature, leading to a damping in the CMB TT power spectrum. 
Thomson scattering also causes additional polarization at scales ~ the horizon at reioniza- 
tion, resulting in excess power at low multipoles in the EE power spectrum. Thus the large 
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angle (small /) EE power spectrum provides valuable information that is complementary to 
the information obtained from the TT power spectrum. 




FIG. 9: (a) shows the large angle EE power spectrum for m x = oo (solid, red), as well as the 
cases with dark matter annihilation (dotted green and blue). Also shown is the first data point 
from the binned BICEP data release, (b) shows the result of 10 5 Monte Carlo simulations for 
multipoles 20 < I < 50, assuming a ACDM model with no dark matter annihilation. The error 
bars are set to n x cosmic variance. The two vertical lines indicate the median value of Xdm f° r 
the m x = 10 GeV model, for n = 4 and n = 3. Only 4.4 (0.25)% of the simulations (with no dark 
matter annihilation) result in a x 2 larger than the median value of Xdm' f° r n = 4 (^)- 

Fig. [9] shows the EE power spectrum plotted for multipoles 20 < / < 50 for the standard 
ACDM model, as well as models which include dark matter annihilation for m x = 10 GeV 
and m x = 1 GeV. Also shown is the first data point from the BICEP results (The WMAP 
experiment provides unbinned power spectra for 20 < / < 50, but the data is prohibitively 
noisy). We have not included data for I < 20 since the power spectrum on those scales is 
significantly affected by standard reionization by luminous sources. For I > 100, the power 
spectrum is damped for the dark matter models, similar to the damping seen in the TT and 
TE power spectra. The current data is insufficient for the large angle EE power spectrum 
to be a useful probe of dark matter annihilation. We therefore perform a number of Monte 
Carlo simulations to quantify the importance of the large angle polarization power spectrum. 

We perform 10 5 Monte Carlo simulations to predict the outcome of a real experiment. 
We restrict our discussion to light dark matter particles of mass m x ~ 10 GeV which are 
motivated by the results of direct detection experiments. We also assume that the correct 
theory is the standard ACDM model without WIMP dark matter annihilation, and test the 
ability of future experiments to constrain dark matter with m x ~ 10 GeV. The assumed 
true parameters are determined by the cosmological model that provides the best fit to the 
TT power spectrum data, and are given by {h = 0.69, n s = 0.97, 10 9 t4 s = 2.245, Q^h 2 = 
0.0225, Q m h 2 = 0.140}. One-step reionization at = 10.5 is assumed. For each value of I, 
the simulated EE power spectrum Cf E is a Gaussian distributed random number with a 
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mean value given by the assumed theory, and a variance 5C EE equal ton x cosmic variance 
(assuming full sky coverage): 

#^ x V27TT ( 21 ) 

We fit each Monte Carlo simulated data set from I = 20 to I = 50 with the assumed 
correct theory. In each case, we compute the value of x 2 (with 31 degrees of freedom). Fig. 
[9^b) shows the number of Monte Carlo simulations that result in a given value of x 2 within 
a bin of size A% 2 = 1, along with the best fit x 2 distribution. We then fit each Monte 
Carlo simulation with the model that includes dark matter annihilation for m x = 10 GeV, 
assuming error bars of 4, and 3 x cosmic variance. The two solid lines show the median 
value of x 2 ; over 10 5 simulations. For n = 4, only 4.4% of the Monte Carlo simulations 
result in a value of x 2 exceeding the median value of Xdm = 45.6 obtained for the m x =10 
GeV scenario. For n = 3, only 0.25% of the simulations result in a x 2 exceeding the median 
Xdm = 57.6. We therefore expect that a dark matter mass m x ~ 10 GeV may be excluded at 
the 95.6% (99.7%) confidence level provided the error bars are smaller than 4 (3) x cosmic 
variance. For comparison, the current WMAP unbinned EE data has an error bar ~ 47 x 
cosmic variance at / = 40. 

It is important to note that the fit cannot be improved by varying A s , the parameter most 
degenerate with the effect of dark matter annihilation. This is because the value of A s is fixed 
independently by the TT power spectrum data. For the case of the standard theory without 
dark matter annihilation (m x ~ oo), the fit to the TT power spectrum data yields Xmm = 
46.8/48 d.o.f. for 10 9 A, = 2.245, with the other parameters set to the values mentioned 
earlier. When the particle mass is reduced to m x = 10 GeV with the value of A s fixed, one 
obtains x 2 — 254/48 d.o.f. which is conclusively ruled out by the data. However, the value 
of A s is not fixed, and increasing 10 9 A S to 2.370 reduces the value of % 2 to 50.8/48 d.o.f., 
which while still disfavored by the data, is not excluded at high significance. Increasing the 
value of A s to better fit the TT data would worsen the fit to the EE data since a damping 
in the small angle TT power spectrum is accompanied by a boost in the large angle EE 
power spectrum. With the new value of A s , the median x 2 for the n = 4 (3) cases increases 
from 45.6 (57.6) to 51.3 (67.8) with 31 d.o.f. 

Let us now consider the Planck experiment. The error bar at multipole / is given by 



SC EE 



Cf E V 2/ + 1 f V 2 

J sky 



1 + (/skyW) 1 m +1)o 2 



(22) 



w = (Tp ix #f whm is the inverse weight per solid angle. cr p i x is the pixel noise, #f w h m is the 
beam full width at half maximum, / S k y is the fraction of sky covered, and Ob = #f w hm/ a/8 ln2. 
The number of pixels w 47r/# 2 whm . The inverse weight per solid angle is then given by 
w~ l = 4ir(AT) 2 / (t x rtboi), where AT is the noise equivalent temperature per bolometer, t is 
the observation time, and nboi is the number of bolometers for the given frequency channel. 
For the 143 GHz polarization sensitive channel of Planck (143P), AT = 82 fiKy/s, n^ Q \ = 
8, and #f w h m = 7' |74J, giving us w' 1 = 2.7(1.1) x 10 _4 /iK 2 f° r ^ = 15 months (3 years). 
Assuming / s k y ~ 0.65 [75], we find that the Planck mission can exclude m x = 10 GeV with 
EE power spectrum data from 20 < / < 50 at ^ 2a with 15 months observation time and 
at > 3cr significance with 3 years observation time. The combined TT + TE + EE data set 
from Planck will provide even better constraints. 
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VI. CONCLUSIONS 



In this article, we have examined how current CMB data can set limits on WIMP dark 
matter annihilation, for the simple models in which the WIMP is all of the dark matter, 
with s— wave dominated annihilation at the thermal rate (<7 a i>) = 3 x 1CT 26 cm 3 /s. Unless 
the dark matter annihilates primarily into neutrinos, one may probe dark matter masses 
m x < 10 GeV using current CMB data. Future data from the Planck mission is expected 
to substantially improve this bound. 

In Section II, we discussed the physics of dark matter annihilation, and computed the 
electromagnetic fraction f em for the annihilation channels XX ~^ bb,cc,T + r~ . We stud- 
ied dark matter annihilation at high redshifts, and obtained an expression for the energy 
absorbed per gas atom, and the ionization and temperature evolution with redshift. 

In Section III, we discussed how the CMB power spectra are modified by dark matter 
annihilation at high redshifts. We performed a maximum likelihood analysis using the 
publicly available CMB Boltzmann code CLASS, and CMB data from the WMAP, SPT, 
BICEP, and QUaD experiments. We obtained the likelihood as a function of WIMP mass m x 
by marginalizing over the cosmological parameters A s , n s , h, Q^h 2 , Q m h 2 . For the XX ~~ ^ bb 
channel, we found that WIMP masses m x < 7.6 GeV are excluded at the 95% confidence 
level, for the simplest dark matter models. We thus find that direct detection experiments 
that prefer a mass m x ~ 10 GeV are consistent with CMB data. We then investigated 
whether constraints on <7g from galaxy cluster observations may improve the bound on m x by 
restricting the range over which A s may vary. Unfortunately, the large degeneracy between 
og and Q m does not allow a precise determination of a 8 , and hence the CMB bound on m x 
is only marginally improved. The bound on m x is strengthened slightly if N e g is accurately 
determined to be larger than the standard model value of N e e = 3.04. For N e g = 3.85, we 
obtain m x > 8.6 GeV (95% confidence, bb channel) which is better than the bound m x > 7.6 
GeV obtained for N cG = 3.04. 

In Section IV, we studied the effect of dark matter halos on the reionization history of the 
Universe. We showed that for optimistic halo parameters, it is possible for dark matter halos 
to substantially reionize the Universe. We obtained an expression for the energy absorbed 
by gas atoms at a redshift z due to particle annihilation at redshift z', and solved for the 
evolution of the ionization fraction as a function of redshift. We then discussed the Ostriker- 
Vishniac (OV) effect as a possible probe of dark matter annihilation at intermediate redshifts 
10 < z < 60. We computed the OV power spectrum for different ionization histories, but 
found the effect to be small except in the most extreme cases. 

In Section V, we examined the importance of the large angle polarization power spec- 
trum as a probe of dark matter annihilation that provides information complementary to 
what is obtained from the temperature power spectrum. Since current experiments do not 
provide high quality polarization data for 20 < / < 50, we performed a number of Monte 
Carlo simulations and showed that the EE power spectrum can help constrain dark matter 
properties. We expect that m x = 10 GeV may be excluded at the 2 (3) a level using EE 
data from 20 < I < 50 provided the error bars are smaller than 4 (3) x the cosmic variance 
value. We expect the upcoming Planck results to provide significant improvements on the 
minimum allowed WIMP dark matter mass. The improved mass bound will be relevant to 
direct detection experiments. 
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